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The scalar-isoscalar propagator of the effective linear cr-model of meson dynamics is investigated 
with the help of an expansion in the number of the Goldstone-bosons. A generic scenario is suggested 
for the temperature/density driven evolution of its pole in the second Riemann sheet. An extended 
temperature range, correlated with characteristic pole locations, is found where the phenomenon of 
threshold enhancement takes place in the corresponding spectral function. 
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I. INTRODUCTION 

O 

The aim of the present investigation is to relate the so-called cr-pole in the second Riemann sheet in the complex 
frequency plane and the behavior of the corresponding spectral function at finite temperature/density. Both objects 
will be determined from the scalar-isoscalar propagator. The calculation suggests a generic trajectory for this pole 
irrespective of the nature of the thermodynamical driving force. We investigated the question of when a well-defined 

, resonance characterized by a Lorentzian shaped spectral function is present. The study is performed in the framework 
of the linear cr-model, used as an effective field theory describing the fluctuations of the chiral order parameter. We 
find the most convenient the application of an expansion in the number of the Goldstone mesons, which is a kind 
of large N approach to the physical excitation spectra of the relativistic 0{N) field theory in its broken symmetry 

' phase. 

The application of the O(N) symmetric scalar field theory to the thermal and finite baryonic density behavior of 
the pion-sigma system was suggested and has been actively pursued for about 15 years, in particular by Hatsuda, 
Kunihiro, and collaborators (for the latest review see Q]). The main physical effect proposed for the scalar-isoscalar 
spectral function is its gradual enhancement near that value of the temperature/baryon density where the phase space 
\ available for the a — > 2n decay is squeezed to zero. 

For the theoretical consolidation of this effect, Hatsuda et al. put forward apparently model independent arguments 
O for the behavior of both the real and imaginary parts of the a propagator. Using an improved version of the finite 
temperature/density perturbative evaluation of the a self-energy O the pole describing the a- resonance moves from its 
' zero temperature/density location smoothly to the location of the two-pion threshold. Both real and imaginary parts 
^ , of the pole location diminish monotonically. The real part approaches the two-pion threshold faster than the imaginary 
k>( ' part vanishes, and the spectral function becomes proportional in this temperature range to the inverse of the imaginary 
part ||. The maximum of the proposed threshold enhancement occurs for that well-defined temperature/density value 
jj^ . when the cr-pole reaches the threshold. 

The effects of partial symmetry restoration realized by the diminishing of /„• was studied also on the unitarized 
7T — 7r scattering amplitude [Q . This quantity was computed and analyzed recently with dispersive techniques in the 
framework of the chiral perturbation theory [^| . 

The approximations made by these authors may be verified soon when reliable spectral functions will be obtained 
using the non-perturbative lattice field theoretical approach || Q|. In the meantime different semi-analytic approxi- 
mation schemes also might shed light on the generality of the proposed arguments. A possible scheme can be based 
on an expansion in the inverse number of the Goldstone bosons. 

The large N expansion has been applied already some 30 years ago to the characterization of the elementary 
excitations of critical O(N) symmetric lattice systems 0|. It was realized that it provides faithful information on 
the excitation spectra of these systems in the full broken symmetry phase between zero and the critical temperature. 

In a recent letter we proposed its application to the present relativistic system, since it avoids all problems of 
principles showing up in other perturbative approaches. First of all its validity does not depend on the rather strong 
self-coupling of the effective cr-model. Second, its result is not sensitive to the choice of the normalization point, 
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that is its change can be always compensated by an appropriate change in the value of the couplings as required by 
the renormalization group. Third, it leads automatically to the same spectra for the elementary cr-field and for the 
scalar-isoscalar quadratic composite field. In other approaches this feature is usually missing or found to be true after 
non-trivial manipulations. Finally, in the chiral limit it provides a correct critical description of the chiral symmetry 
restoration. 

The leading order large N approximation has been applied to the Goldstone boson scattering by Chivukula and 
Golden p~T| ] at zero temperature. They have explicitly checked that this approximation fulfills in the scalar-isoscalar 
channel the unitarity condition and also satisfies the Adler-zero condition. 

The large N leading order amplitude (N — l)A(s) has been compared also with the existing phase shift data in 
the I = J = channel of the ir — n scattering by Dobado and Morales |Q . They have completed the leading order 
amplitude by subleading terms, dictated by the requirement of crossing symmetry: A 00 (s, t, u) — 3A(s) +A(t) + A(u). 
Though this is not a systematic next-to-leading order computation the authors found a satisfactory fit to the relevant 
phase shift <5§(s) upto^/iw 600MeV. 

To our best knowledge no application of the large N approach at finite temperature/baryon density was attempted 
to date. In view of the importance of t- and u-channel exchange contributions to the ir — ir scattering amplitude, 
emphasized in the recent literature Jl3| , a fully satisfactory large N treatment will also require the calculation of the 
next-to-leading order contribution to the a propagator. The leading order calculation presented in this paper is a 
necessary intermediate step also towards this goal. 

The fact that for T — the real and imaginary parts of the er-resonance are of comparable magnitude implies that for 
its finite temperature description we cannot restrict ourselves to the immediate vicinity of the real axis. An exception 
is when the er-resonance gets close to the two-pion threshold, which is the situation where threshold enhancement 
develops. For the explanation of the detailed features of this phenomenon it is unavoidable to trace the complete 
temperature driven pole trajectory. The main purpose of the present paper is to carry out such an investigation 
within the leading order of the large N expansion. According to our calculation to be presented in this paper, the a 
self-energy continued analytically into the lower halfplane leads to a pole trajectory (when the temperature T or the 
baryonic density, ub is varied) whose real part assumes smaller values than twice the pion mass, while the decrease 
of the imaginary part is not drastic. As a consequence of this the pole looses its meaning as a resonance and at 
the same time in the spectral function the direct neighborhood of the threshold will be emphasized. When further 
increasing T/ub the trajectory hits the real axis on the second Riemann sheet and moves along it before it reaches the 
threshold. This pole evolution will be demonstrated to be generic in the sense that it is insensitive to the variation of 
the parameters of the theory. One can argue that it can stay valid beyond the leading order large N approximation, 
since it is unlikely that the trajectory would hit the threshold point directly. Only in the limiting chiral symmetric 
case do we find the smooth behavior which was proposed by Hatsuda et al. . 

One should be aware of the fact that in scalar models a tachyonic pole is always present |[4| |l6| related to the 
Landau-ghost phenomenon. It restricts the range of variation of the renormalizcd parameters where the model can be 
used in an effective sense at all. For this reason we find in the leading large N approximation that a parametrization 
accounting for the T = phenomenological data of the ir and a mesons can be achieved only approximately. However, 
we shall argue that the relationship between the pole trajectory and the variation of the spectral function to be 
described below might not change qualitatively. 

The presentation of the paper is aimed at a self-contained, technically transparent description. In section II the 
Schwinger-Dyson equations for the finite temperature two-point functions of the linear a model are given to leading 
order in N . Their analytical continuation onto the second Riemann sheet in the complex frequency variable is 
presented in section III. Its explicit expression was found by studying the bubble diagram, describing the splitting of 
the cr-field into two off-shell pions. In section IV we analyze the temperature driven evolution of the a pole-trajectory. 
This investigation makes use of the T = parametrization of the model as an input, therefore one first works out the 
details of the physical and unphysical poles for T = 0. Here we fix all couplings in a way to achieve the closest possible 
characteristics of the a to the particle data, and still staying by a factor of 2 to 3 below the "energy scale" of the 
tachyon. In the second part of the section we shall argue that the pole-trajectory found by varying the temperature 
remains qualitatively the same when the system evolves under the variation of the baryonic density. A detailed 
discussion of the change in the pole trajectory pattern with the pion mass is also presented. The spectral function 
of cr is computed in section V by approaching the real axis from the physical upper half-plane. We shall analyze 
also the function which arises when the threshold factor is divided out. It will be demonstrated that in an extended 
temperature range the maximum of the spectral function is located in the closest neighborhood of the (T-dependent) 
position of the two-pion threshold. The extension of this interval is very well correlated with the piece of the pole 
trajectory, when its real part is below the two-pion threshold. In section VI the conclusions of the present study are 
summarized. 
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II. LEADING LARGE TV EXPRESSION OF THE PROPAGATORS AT FINITE TEMPERATURE 



The appropriate parametrization of the O(N) symmetric Lagrangian for a large TV expansion has the following 
form: 



[d^ a d^4> a - m 2 (j) a (j) a } 



A 



24TV V 



AQ ) 2 (/) 2 + ViV/ l 1 . 



(1) 



The last term explicitly breaks the O(N) symmetry and introduces non-zero mass for pions. 

In the broken symmetry phase one separates the expectation value $(T) of the field, which points along the direction 
a = 1 in the internal space 

(2) 



In the following all quantities will be computed to leading order in the large N limit. 

The quantum fluctuations of the order parameter are divided into a longitudinal mode, which represents the a 
meson and the transversal ones. The latter correspond to the Goldstone modes, the pions. Their mass, mc{T) is 
determined as the pole of the resummed pion propagator, in which the tadpole contribution is calculated with the 
pion propagator selfconsistently. This results in the following gap equation: 



A 



A 



mUT) = ™ 2 + -^(t) + — dry) 



A 



$ 2 (T) 



cPk 1 
(2tt) 3 2^ 



(1 + 2n(w*)) 



(3) 



where n(wfe) = l/(exp(a;fe/T) — 1) and Uk = \Jk 2 + mg(T). 



The equation of state obtained from the requirement 



is as follows: 
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m 2 G 



(T) = h/^(T). This is 



Comparing this with Eq.(|j|) one can observe that for ^ consistency requires 

precisely the Goldstone theorem in the presence of explicit symmetry breaking. 

At leading large TV order the contribution to the longitudinal self energy is given by the sum of the contributions 
corresponding to the diagrams of Fig.[j], i.e. the bubble series. On both internal lines of a bubble exclusively the 
propagation of pion fields are taken into account to leading order in TV. In each term of the bubble contribution 
there is a common multiplicative vertex contribution coming from the two edges of each diagram. The bubble series 
consistently takes into account the non-zero classical value for ^(T), partly by the implicit dependence of m G (T), due 
to the gap equation and also by the two legs of the effective 4-point vertex formed by the sum of the bubble-series. 




oo-o- 



FIG. 1: Diagrams determining the self-energy of the a field represented to leading order in TV. The external solid lines 
correspond to the a, while in the internal bubbles pions propagate. The dashed line represents the expectation value $ = $(T). 
The vertices can be read from Eq. {jlh after the shift defined in Eq. (kh is performed. 



By adding the tree level mass m 2 to the self-energy of the a field determined in the background $(T), using Eq. (^) 
one finds for the effective cr-mass a simple expression in terms of b{p) , denoting the value of the single bubble diagram: 
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The bubble contribution with external momentum p = (po, p) is the sum of a zero temperature and a T-dependent 
part, b > (p) — (p) + b^(po,p). The superscript > hints at the fact that the expression of the bubble contribution is 
valid in the upper p half-plane. The explicit expressions of the two terms read as follows: 
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where n, = l/(exp(/3o; i ) - 1) and = (q 2 + m^T)) 1 / 2 , u 2 = ((q + p) 2 + m 2 ,^)) 1 / 2 , and e > 0. 



Us 



ling cut-off regularization the zero temperature bubble contribution is 
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It is common to discuss the spectral function for p = 0, when the expression of the finite temperature bubble 
contribution simplifies to: 
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In view of the quadratic and logarithmic cut-off dependencies which appear in Eqs.(|^),(^) and (||) a mass- and 
coupling constant renormalization is necessary. It requires the introduction of a normalization scale M a . This scale 
should lie below the scale of the tachyonic pole to be discussed below. Its choice within the relevant range where 
the O(N) model can serve for the effective description of hadron dynamics should not affect sensitively the physical 
results. 

The expressions of the renormalized couplings go beyond the accuracy of the usual one- loop relations. Actually, 
they ensure that a change in the normalization scale can be compensated in all formulae below by an appropriate 
change in the couplings. The following non-perturbative mass- and self-coupling renormalizations are introduced: 
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In terms of the renormalized couplings the equation of state Eq.(^|) can be cast into the following explicitly finite 
form: 
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ma(T)/T exp(y) 

value of the Goldstone mass and of the expectation value of the a field, 



Here 777^0 an d ^0 stand for the T 
respectively. 

The formal expression of the effective mass term of the cr-propagator is unchanged after renormalization, just 
Xr should replace A and for the T — contribution of b(p) in Eq. (||) the scale M5 2 is put in place of eA 2 . 

The finite part of the zero temperature bubble contribution defined in this way has different forms depending on 
the range of po values: 
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III. ANALYTICAL CONTINUATION OF THE PROPAGATORS ONTO THE SECOND RIEMANN 

SHEET 

The interpretation of the temperature driven variation of the spectral function will be based on the study of the 
variation of the scalar-isoscalar pole mass, that is the zero of G~ 1 (po, p = 0) = p\ — m 2 in the lower po halfplane. For 
this it is necessary to construct an analytical continuation of the longitudinal propagator onto the second Riemann 
sheet. By Eq. (j5j) it is clear that the problem is equivalent to the continuation of b{po), which will be discussed in this 
section. 

We decided to perform the analytical continuation in such a way that b(po) varies continuously when the real axis is 
crossed above the two-pion threshold po > 2ma{T). This implies that the bubble contribution is discontinuous across 
the real axis for —2mc{T) < po < 2tt7g'(T). We will actually see that there it is Re b(po) which is discontinuous. 

Above the threshold on the real axis both the real and imaginary parts of the zero temperature bubble are continuous 
as one can see from the second line of Eq. (O) . In view of this, for values of po below the real axis (lying on the second 
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Riemann sheet) we simply use the expression written in the second line of Eq. (|12|). We denote this continuation 
6q (po). (In this sense the indices ">" and "<" on 6 are redundant, yet we keep them for the sake of clarity.) 

We turn now to the finite temperature part of the bubble. For real values of po, above the threshold the real part 
of 6t(po) can be obtained by taking the principal value in the right hand side of Eq. (|^): 



i „ r ,yx 2 -m' a m/Ti 1 



te *w-g p ^« ^ii7aeg expM _ r (13) 

Introducing the parametrization p Q = Rep + is we can evaluate the imaginary part of b^{po) m the neighborhood 
of the real axis both for e > (the physical prescription) and for e < 0: 
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With this one obtains: 
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In order to ensure the continuity of the imaginary parts an extra term has to be added to the expression used in 
the upper half-plane for the bubble: 



i , , „ / 4m%(T) 

to"™** f— ft 



b<( Po ) = b>( Po )-i-n(p /2)\ll <§-L. (16) 



For later use (among others for the computation of the spectral function) it is useful to write explicit expressions 
for the physical values of b(po) on the real positive axis. We have b^,(po) — Re b^ipo) + *I m b^(po), where the first 
term is given by Eq.(^), and 



Im b>( P0 ) = -—* 1 ^n{p /2)Q (p - 2m G {T)) . (17) 
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The bubble is fully real below the threshold and it is clear that the integral in Eq. ([l3]) is not singular for po < 2m G (T). 
Its integrand has for po = 2mg(T) an integrable square root singularity: 

1 f°° 1 1 
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For complex po also in the lower halfplane one can use for bg(po) the expression written on the second line of Eq. 
(O) and the complete expression of b^ipo) reads as 
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IV. THE TEMPERATURE DEPENDENCE OF THE cr-POLE 



In this section the temperature driven variation of the pole of G a (po), located in the fourth quadrant of the complex 
po plane will be found. We shall argue that the generic scenario of its variation is realized independently of what 
thermodynamical quantity would drive this variation. As an illustration of this we shall discuss the pole trajectory 
under the change of the baryonic charge density implemented following Hatsuda, Kunihiro, and Shimizu jj). 

For the solution of the equation 

G-\ P o) (l ^&<(po)) = (p§ - ml{T)) (l - ^(b<(p ) + 6<(p ))) - ^$ 2 (T) = (19) 

one obviously should know <3>(T) (and m G (T) — h/Q(T)), calculable from the renormalized equation of state jll]). 
This equation requires two inputs. The phenomenological (physical) input is $§/mg. = $o//i = /^/4mg ~ 0.11, 
but in the equation of state also the normalization scale appears explicitly. 

In order to simplify the formulae we have decided to choose the absolute value of the pole location at T — for the 
normalization scale Mo- Therefore one has to find first the er-pole for this temperature and only then one can turn to 
the discussion of the finite temperature variation. 
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A. The poles of G CT (p ) at T = 



We parametrize the solution of the equation 

g-Hpo) (i-^Cpo)) 



(20) 



in the form po — Moexp{—iipo), < ipo < ir/2. Note, that the renormalization scale is fixed in proportion to Togo 
once a renormalized coupling is chosen. 

One can introduce instead of M and ipo a more convenient parametrization pq = 2mQ + Mq exp(— i(f>o), where Mq 
and <po are uniquely determined by Mq and tp^. Attention has to be payed when switching from one parametrization 
to the other because changing the scale Mq means changing the renormalized coupling constant Xr. 

In addition one has to deal with care when realizing numerically the analytic continuation. For the square root and 
for the evaluation of the argument of the complex numbers appearing in 6^ (p) one has to choose a phase convention 
which ensures the continuous variation of the complex phase of the final complex number with ipo . Good guidance 
for how to define the argument of complex numbers in the process of evaluating complicated multivalued functions is 
provided by their series expansion near the positive real axis above the threshold. 
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FIG. 2: The real and imaginary parts of the physical poles at T = 0. Also shown is the logarithm of the tachyon pole position 
in proportion to the mass of a for various temperatures. The lines appear in the same order downward from above on the right 
side of the figure as the labels in the key. 



The output is Mq/^q and ipo in terms of which we obtain the mass M a = Mocos</5o and the width T = 2Mo8W.(po 
of the a as the real and imaginary parts of the pole. These are shown as a function of the renormalized coupling 
Xr in Fig.||. At Xr = 400 one finds the ratio M a /T ~ 1.4 with M a = 3.95/^- These values are away from the phe- 
nomenological expectations [[l7| [ll| , but these are the best values we can reach in the leading large N approximation. 
The use of higher values of Xr which appear to be closer to the observed numbers might not be advisable since the 
a tachyonic pole described below comes very close to the scale Mq for that coupling region. We will use Xr = 400 
henceforth in the finite temperature calculations. 

Scalar theories are known to have a tachyonic pole related to the Landau-ghost that is a zero of the inverse 
propagator on the positive imaginary axis po = iM^. Equation ( ^p|) takes the following form when looking for M^: 
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With a suitable parametrization, for example, Ml /Mq 
$o/ m Go an d using the value of M /$o obtained by solvint 
restricts the range of validity of the theory is shown in Fig 
non-zero temperatures. 



= exp(z), one can solve Eq.@ for z at a given value of 
Eq. (|2fj) . The logarithm of the ratio Ml jM a whose value 
3 as a function of Xr, not only for T — but also for some 
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B. Finite T/ns behavior of the a-pole 

The features of the numerical solution of Eg. (|l9|) will be discussed in the main part of the present subsection. It will 
be pointed out that when moving on the second Riemann sheet its root approaches and eventually hits at a certain 
temperature the unphysical real axis below the two-pion threshold. For the determination of the poles on this piece of 
the real axis one can evaluate the analytic functions directly. The real solution provides a useful check of the solution 
based on the complex equation. We give here the corresponding formulae explicitly. 

In both parametrizations of the pole one has on the real axis below the threshold (po < 2mn(T)), w. — ir, (p = it, 
respectively. Then one has (1 — 4toq(T)/pq) 1 / 2 = — i(4mg,(T)/pg — l) 1 / 2 = — iQ. Using this in Eq. (12) one obtains: 



4m%(T) 1- y/l - 4mUT)/p 2 , . \ _ , n A Q^ 1 + * 



-i i in ===== ,„ , = iQ In -— — — 

Pi \ 1 + v/1 - 4m£(T)/pg J V Q 

= -2Qarccot(g _1 ) - irQ = 2Q(arccot(Q) + tt). (23) 
So, on the real axis of the second Riemann sheet, below the threshold one has: 
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(Here we use a somewhat redundant superscript "~" on b < (p) in order to emphasize that the formulae refer to a 
continuation below the threshold). 

The evaluation of b > (p ) below threshold on the real axis becomes relevant for temperatures when the pole "climbs 
up" from the second onto the first Riemann sheet and represents a stable a particle. When approaching this portion 
of the real axis from the physical upper halfplane one trivially finds the corresponding real expression. For real po the 
equations which determine the poles on the physical and unphysical real axes below the threshold (po < 2rriG(T)) are 
the following: 



(Pi 



1{T)) (l - ^ (b>( P0 ) + 6>( P0 ))) - ^$ 2 (T) = (26) 
and 

(pi m%(T)) (b<~(Po) + &£-(Po))) - ^ 2 (T) = 0, (27) 

respectively. 

In Fig.|| the trajectory of the real part of the pole position is shown as it evolves with the temperature. In the 
same figure also the T-dependent location of the two-pion threshold appears. One notes that Re(po) crosses below 
the actual position of the threshold for T** /mQQ ~ 0.68. As it is seen from Fig.|]lm(po) hardly diminishes until this 
temperature is reached, therefore the relative broadening of a actually increases. Similar conclusions where drawn in 
the context of chiral perturbation theory in |19[ . Above T** the imaginary part decreases faster and the pole position 
is landing on the real axis for T rea i / mco ~ 0.93. 

One can study the solutions of Eq.(|l9|) also in the upper po halfplane on the second Riemann sheet. One finds that 
starting from T — 0, there exists a "mirror" -root, which arrives at the same point of the real axis for T = T rea i. This 
collision of the poles results in two oppositely moving real solutions for higher temperatures. The solutions of Eq. ( p7| ) 
fully confirm this scenario. 

The pole moving upwards catches up with the threshold for T*/rriGQ — 1-074. (The other pole first moves down- 
wards, later its motion changes direction, but it lags behind the position of the threshold in the whole temperature 
range of interest, see. Fig. |^.) It does not stop there, but moves further with increasing temperature, now on the real 
axis of the physical Riemann sheet. This part of the scenario is confirmed also by the direct solution of Eqs.(|27j) and 
(p6|). The T-dependent position of the stable physical a particle is also displayed in Fig.||. One has to notice that all 
scales increase with the temperature, therefore the tachyon pole puts a strict temperature limit to the validity of the 
proposed effective treatment of the pion-sigma system. 

The scenario obtained is clearly different from the one suggested by Hatsuda, Kunihiro, and their collaborators [0. 
In their various approximate descriptions the real part of the pole position never goes below the two-pion threshold 
until it has a finite imaginary part. In order to test the generic nature of the pole evolution found above, we have 
followed the procedure proposed by Hatsuda et al. B for the introduction of finite baryonic charge density tib into the 
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FIG. 3: The temperature dependence of the real and imaginary parts of the cr-pole 



effective pion-sigma dynamics. In summary, the nonzero value of ng results in a rescaling of the vacuum expectation 
value $ P0|| - At T = 0, for low densities they propose 

*(n B ) = (1 - Cn B )$ , (28) 

with C = 0.2 — 0.3. With the parametrization po = r]Moexp(—iip) we have solved Eq.(|20|) for rj and (p with $(ns) 
replacing <!>n everywhere. This equation is much simpler, since the details of its continuation onto the second Riemann 
sheet are self-evident. Still, we find qualitatively the same pattern for the pole trajectory as one can see in Fig.|. In 
another test the strength of the explicit symmetry breaking (h) was gradually decreased. The distance of the point 
where the complex solution arrives to the real axis from the two-pion threshold monotonically decreases with the 
decrease of h. In the chiral limit it approaches smoothly the origin as it should for a true phase transition. It is 
worthwhile to point out that the most recent dispersive investigation of Yokokawa et al. || investigates exactly this 
limit. They find the same smooth behavior, therefore there is no conflict between the results of the two approaches 
yet. It will be interesting to see the effect of explicit symmetry breaking in their approach. 
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FIG. 4: The dependence of the real part of the pole position on the variation of due to non-zero baryonic density 

Finally, we have studied systematically the deformation of the T-driven pole trajectory when 2mco/ fir is decreased 
gradually with f n kept constant. A quite interesting pattern appears in Figj^. The trajectory reaches closer to the 
negative imaginary axis as the strength of the explicit symmetry breaking diminishes, before it turns to the real axis, 
and eventually ends at the two-pion threshold. At some value 1.317 > 2mGo,imag/ fir > 1-316 it touches first this axis, 
but is "reflected" from it back into the fourth quarter. For higher temperatures it will have again a non-zero real 
part. 

In order to understand what happens, it is convenient to search directly for poles on the negative imaginary axis. 
It turns out that already for 2togo/ fir — 3.01 one has an infinite number of approximately equidistantly located poles 
along this axis for finite temperature. Below we call the pole located the closest to the origin the er*-pole. Already 
at low temperature the distance of all these poles from the threshold is much larger, than that of the cr-pole. This 
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FIG. 5: Trajectory of the complex cr-pole for various values of 2mao/ fw Note the tendency of the trajectory to approach closer 
to the imaginary axis as moo decreases. 

explains why the latter dominates near T = the behavior of the spectral function, as will be argued in the next 
section. The distance further increases with the increase of the temperature. 

It turns out that for the above quoted value of 2rriGo,imag/ f-w the cr-pole and its mirror from the third quarter touch 
the imaginary axis exactly at the location of the highest (negative) imaginary <r*-pole at that temperature. The result 
of the "cr — cr*" collision is the reflection of the two complex poles back into their respective quarters. 

Further decreasing 277iGo// 7r the pair of mirror poles arrives onto the negative imaginary axis below the highest 
imaginary pole. (In Fig. o it is above cr*, since the negative imaginary axis is directed upwards.) The colliding 
poles now give rise to a purely imaginary pair, one member of which moves towards the origin the other one moving 
the opposite direction. The pole moving towards the origin collides at some higher temperature with the oppositely 
moving genuinely imaginary cr*-pole, and they are pushed back into the complex quarters as mirror poles. Eventually, 
the pole in the fourth quarter will land on the real axis and moves up to the threshold for T*(togo)j where it is 
converted into a stable particle pole on the physical sheet, as described earlier in this section. In view of the multiple 
pole collisions it is clear that the stable high-T a particle is not directly related to the T = complex cr-pole. 

In the chiral limit the cr-pole reaches the origin and describes the phase transition restoring the chiral symmetry 
of the model. This simple trajectory has been already discussed in our previous publication JlC[ |. The presence of 
explicit symmetry breaking with realistic strength led to a rather spectacular change in this scenario. It is a valid 
question whether a smooth continuous deformation of the pole trajectory connects the case of the explicit symmetry 
breaking with the chiral limit. 

In the region, when the pion mass is much smaller than the temperature for which the complex pole becomes 
purely imaginary, one can find analytically the first cr*-pole in the infinite sequence described above. Its location 
is given approximately as p 2 ,, ~ m^(T)T / <& 2 (T) . In view of the fact that T/$(T) ~ 0(1), this pole goes faster to 
zero than the pion mass. Our numerical study shows that the point where the pole lands on the negative imaginary 
axis does not change more than 10% between the pion mass rriGo,imag and the chiral limit. Therefore the interval 
of temperatures for which the pole moves on the imaginary axis increases with decreasing pion mass (see Figj^). 
Eventually for mco = the highest negative imaginary pole stays (with zero residuum) in the origin and the scenario 
characterizing the chiral limit sets in smoothly. 

For symmetry breaking much smaller than the parameter characterizing the onset of the dynamical scaling in the 
chirally symmetric case one can even experience the realization of the scaling behavior. 

To some extent the above complicated trajectory is to be expected, since it is highly "improbable" that the roots of 
a complex equation would move smoothly to a specific point (e.g., the two-pion threshold) of the real axis, irrespective 
of the variation of its parameters. The limitations of the true resonance interpretation of the cr-pole will be discussed 
in the next section, when the correlation of its location with the measurable spectral function p a will be discussed. 

V. THE T DEPENDENCE OF THE SPECTRAL FUNCTION p a 

The spectral function of the order parameter field a is defined using the expression of the propagator in the physical 
half-plane as: 

Pa (p ,p,T) = — lim ImG a ( PQ +ie,p,T). (29) 

7T e— >+0 
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The leading order large N expression of spectral function, at p — 0, is given by 

Pa (po,0,T) = A|$ 2 (T)Im6>(p )/187r x 



ml{T)) (1 - ^Re6>(p )) - Af $ 2 (T)] 2 + (p 2 m G (T)) 2 ^(Imb>(p )f 



(30) 



We have seen that lmb > (po) ^ only for po > 2itig(T), if there is no stable particle pole below the threshold. 
Therefore the spectral function is nonzero only for po values above the threshold until T < T*. For this reason the 
numerator and the second term in the denominator of Eq. (|30|), that is Imb > (po) goes to zero when po — > 2mc(T) + 0. 

The scalar-isoscalar spectral function is displayed in FigM One notices the shift of its maximum towards the 
two-pion threshold with increasing temperature, though its width does not decrease initially. For T > T** the shape 
of p @ becomes cuspier. Finally, a very high value of the maximum is experienced numerically around T/ttiqq ~ 1.07. 
Above this temperature the value of the maximum gradually diminishes and its location shifts increasingly farther 
from the threshold towards larger pq. 

For a qualitative interpretation of p a one verifies that for the temperature T* w 1.074m<3o, the first term of the 
denominator of Eq.(^0|) vanishes at the threshold po = 2mc(T*), that is Eq.(p6|) is equivalent to the condition for the 
vanishing of the first term in the denominator of Eq. (|3C 
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= 0. 



(31) 
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FIG. 6: The scalar-isoscalar spectral function at various temperatures 



Because the term containing the real part of the bubble vanishes more rapidly at T* as po — > 2ma{T*) than the term 
containing the imaginary part of the bubble, the behavior of the spectral function around the threshold is dominated 
by the imaginary part of the bubble Paipa, 0, T*) ~ l/lmb > (p ) = 1/(1 — Am 2 G {T*) jp\) x l 2 . This is formally the same 
behavior exploited by Hatsuda et al. B when arguing in favor of the generic nature of the threshold enhancement 
phenomenon. It is worthwhile to emphasize, however, the obvious fact that according to our calculation the spectral 
function has nothing to do with the imaginary part of the second sheet pole which is a purely real (unphysical) 
quantity in this temperature regime, and the singular behavior is due to the fact that the pole along the unphysical 
real axis moves towards the threshold when T — > T* . 

One might suspect, that the coincidence of the pole position at T = T* with the threshold would lead to a stronger 
~ S(po — 2m G (T*)) singularity for this temperature. This is not true since one can easily demonstrate the vanishing of 
its residuum for T — T* . The residuum of the stable a pole appearing on the physical sheet for T > T* continuously 
increases with the temperature. 

The near threshold enhancement of the spectral function is maximal at T*. When T > T*, the position of the 
maximum moves away from the threshold and its height diminishes. This is consistent with the requirement arising 
from the sum rule, J dpoPoPo-(Po) — 1 m presence of a stable cr-pole with increasing residuum. 

The enhancement sets in gradually and for T > T** the maximum of p a stays very close to the actual threshold 
position. We can argue rather convincingly for a certain physical significance of 7 1 **, when displaying p\{p^) = 
p a (jpo) / yl — 4 : m 2 j(T)/pQ. The argument for this operation is the fact that the phase space volume is just proportional 
to the factor divided out. In Fig.|7j one can follow the position of the maxima of p„ and p\ relative to the two-pion 
threshold as a function of the temperature. One sees that the position of max(p (T ) touches the threshold only in a 
single point T = T*. On the other hand the position of max(pJ.(T)) approaches 2ma(T) rather steeply, but starting 
from T = T** the distance is found numerically always smaller than 10~ 5 . 
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FIG. 7: The temperature dependence of the locations of the maxima of p a and of p\ as compared to the two-pion threshold 

In conclusion of this section one sees that the qualitative changes in the spectral function can be well interpreted 
with the help of the scalar-isoscalar pole located in the lower p halfplane. The threshold enhancement occurs in 
an extended temperature region, T** < T < T* . Any analysis attempting the reconstruction of the in- medium a 
resonance from some experiment is necessarily based on the behavior of the scalar-isoscalar spectral function. From 
the above discussion it is suggestive that by the enhanced signal coming from the neighborhood of the threshold one 
would be led to the conclusion that in the temperature range T** < T < T* the "er" moves together with the two-pion 
threshold. But no trace of any Lorentzian resonance shape can be detected in the spectral function. 

VI. CONCLUSIONS 

In this paper we have presented all the technical details of analyzing the pole trajectory of the propagator de- 
scribing the fluctuations of the chiral order parameter in the linear sigma model with the help of a leading large N 
approximation. It is worth noting that the propagator of the composite field (0 a (x, t)) has the same poles as that of 
<j> in the broken symmetry phase (see Ref. ]l(| and references therein). Detailed mapping of the trajectory under the 
variation of the temperature/baryonic density as well as the parameter controlling the explicit breaking of the chiral 
symmetry leads us to conjecture that generically the real part of the complex pole will decrease below the two-pion 
threshold energy with a finite imaginary part. The application of higher order corrections in the large N expansion 
are expected to give no more than 25% correction. Therefore only quantitative change may happen in the scenario 
presented in our paper. 

Our main result is the prediction that in a whole temperature range T** < T < T* the maximum of the spectral 
function is rather close to the actual location of the two-pion threshold. In this sense one can speak about an interval 
of threshold enhancement. The point of maximal enhancement is T = T* . 

The usual connection between the second Riemann sheet pole and the spectral function is restricted to the tem- 
perature range T < T**. In this region the spectral function has an approximately Lorentzian shape, its maximum is 
centered nearly at the location of the real part of the pole, its width scales with its imaginary part. When approach- 
ing T** from below, the approximately Lorentzian shape of the spectral function will be distorted, in particular it 
is loosing its symmetry. With considerable compromise one can extend the cr-particle interpretation of the second 
Riemann sheet pole up to T** . 

The complete loss of this characteristic means that no resonance interpretation can be given to the pole in the 
T** < T < T* range. On the other hand any phenomenological analysis of a particle signal proportional to the spectral 
function will be peaked near the two-pion threshold in the above temperature interval. The real-time propagator in 
this channel will be dominated by the contribution of frequencies slightly above Ivtlq (T) . In this sense an increasingly 
narrow cr-signal can be detected at u> ~ 2mc(T) in the final two-pion spectra. 

Above T* it is a stable physical cr-particle, which reappears. Accordingly, the peak of the spectral function dimin- 
ishes. By the increased effect of the phase space factor its maximum will be reached again at frequencies farther above 
the two-pion threshold. 
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